arXiv:astro-ph/0512212v3 24 Aug 2006 


Prospects of LIGO for constraining inclination of merging compact binaries associated 
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We study prospects of a method to constrain the inclination of a coalescing compact binary 
by detecting its gravitational waves associated with a three-dimensionally localized (direction and 
distance) short-hard gamma-ray burst. We take advantage of a synergy of these two observations, 
and our method can be applied even with a single interferometer. For a nearly face-on binary the 
inclination angle I can be constrained in the range 1 — SNR -1 < cos 7 < 1 {SNR: the signal 
to noise ratio of gravitational wave detection), provided that the error of the distance estimation 
is negligible. This method would help us to study properties of the short-hard bursts, including 
potentially collimated jet-like structures as indicated by recent observation. 


I. INTRODUCTION 

The gamma-ray bursts have been known to be divided into two classes, the long-soft bursts and the short-hard 
bursts (SHBs). While the former are likely to be produced at explosions of massive stars in star forming galaxies 
typically at high redshift z ^ 1 Q|, the nature of SHBs has been a long-standing mystery. However, recent discoveries 
of X-ray afterglows of SHBs by Swift and HETE satellites allowed us to localize them accurately and rapidly enough 
to specify their host galaxies and finally determine their distances IlSiii 

One of them, GRB 050724 was found in an elliptical galaxy at z = 0.257 with an old stellar population Q, and 
GRB 050509b is likely to be in a similar galaxy at z = 0.225 Q. While GRB 050709 was in a star forming galaxy at 
2 = 0.160, its light curve excluded a supernova association Q. These results support that coalescing compact binaries 
(double neutron stars (NS+NSs) or black hole-neutron star (BH+NS) systems) are the promising origins of SHBs [?|, 
though the estimated typical age of these binaries for SHBs are longer than that of known NS+NSs in our galaxy 

il 

Coalescing compact binaries are also promising sources of gravitational radiation for LIGO and other ground-based 
interferometers [lQ. A one year scientific run (S5) is ongoing with LIGO that has sensitivity to detect NS+NSs to 
~ 15Mpc El- Recent theoretical analysis predicts that the probability of a simultaneous detection of gravitational 
waves by LIGO and a SHB by Swift in one year is ~ 30% for BH-NS merger and ~ 10% for NS+NS (depending 
on the lower end of the luminosity function of SHBs) [1]. Therefore, we might soon experience the first detection of 
gravitational waves associated with a localized SHB. 

The observed afterglows of two SHBs showed steeper power-law decays that indicate SHBs have collimated jet-like 
structures Hi, as found with long-soft bursts. The estimated beaming factor is ~ 0.03 for GRB 050709 |5] and 
~ 0.01 for GRB 050724 [§]. Here we defined the beaming factor as the fraction of 4-7 t steradians into which jets are 
emitted. This enabled us to estimate their total energies ~ 3 x 10 48 erg that is smaller than the long-soft bursts by 
~ 2 orders of magnitude. If a SHB is associated with a coalescing binary, it is likely that the orientation of the jet 
is aligned with the angular momentum of the binary that would be clearly imprinted on the observed gravitational 
waveform. 

In this paper we propose a method with which the inclination of a binary will be interestingly constrained as a 
synergy of three-dimensional localization by electro-magnetic waves (EMWs) and observation of gravitational waves, 
even using a single interferometer. Therefore, LIGO could provide us an important geometrical information to 
understand properties of SHBs. 

This paper is organized as follows: in §11 we describe our method to estimate inclination of binaries with using a 
single gravitational wave interferometer. Expected error for our method is studied in §111. Then, in §IV, we extend 
our study for observation with multiple detectors, such as, LIGO-VIRGO network. Brief discussions are presented in 

§v. 


II. GWS AND CONSTRAINT FOR INCLINATION 

For simplicity we use the restricted post-Newtonian description H2| with neglecting precession induced by spin 
that might be important for BH-NS, but not for NS+NS m (see also Il4l| for recent analysis). The two polarization 
waveforms in the principle polarization coordinate are given as mm 

h+(t) = Bf 2 ^ 3 {l + V 2 ) cos($(i)), 


(1) 
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FIG. 1: Configuration of a binary and an interferometer (X and T-axes; directions of its two arms). The unit vector N 
represents the direction of the binary, and the unit vector L is the orientation of its angular momentum. The angle / is the 
inclination with cos I = L ■ N = V. The polarization angle ip specifies the direction of L around the vector N. 


h x (t) = R/ 2 / 3 (2R)sin(<f>(f)), (2) 

where <b(f) = 27r f* f(t', M c ,..., ti)dtt' + ip c is the phase of the waves, and tj and y>c are constants. The parameter V 
is defined by 

V = cos I (3) 


with the inclination angle I shown in figure 1, and is the primary target in this paper. The chirp mass M c is the 
most important parameter to characterize the time evolution of the frequency / = 1/27T d^/dt and given by two 
(redshifted) masses of the binary as M c = mj /,5 mf 5 (mi+m 2 ) -1 / 5 . The intrinsic amplitude B at the quadruple order 
is given by the chirp mass M c and the distance (more precisely, luminosity distance) r to the binary [l0|, [l(j, 


D o G 5 / 3 M c 5/3 tt 2 / 3 
B = 2 —^— 


( 4 ) 


Response h(t) of the interferometer to the two polarization modes is written in terms of the beam-pattern functions 
F + and F x as 


h(t) = F+h+(t) +F x h x (t). (5) 

The beam-pattern functions are determined by three angles 0, (j> and ip as 0 

F + = ai cos 2 ip — 02 sin 2 ip, (6) 

F x = ai sin 2ip + 02 cos 2ip (7) 

with ai = \ (1 + cos 2 6) cos 2 <p and 02 = cos 9 sin 2 <p. The angles (0, cp) represent the direction of the binary in the polar 

coordinate attached to the interferometer as in figure 1, and ip is the polarization angle that fixes the axial direction 

of the vector L around the direction N (see e.g. figure 1 in (l3j|). 

When we take Fourier transformation of the response h(t), the result is formally given as 

kf) = AR(V, ip, 9 , cP)r 716 exp[itf (/)]. (8) 

Here we used the stationary phase approximation, and the phase 'F(/) is a real function. The overall amplitude of 
the signal AR(V, ip , 9, (p) is given by llTl 


R=\J{(V2 + l)F + } 2 + (2VF x )i , 


( 9 ) 
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TABLE I: Determination of Q = (AR)/AR max 


notation 

observation 

error 

A 

M c \ GW (chirp) 

~ 0 

(oc Mp /0 /r) 

r: EMW (localization) 

AH 0 /H 0 , velocity 

(AR) 

GW (amplitude) 

~ (SNR)- 1 

Rmax 

(9,cp)-. EMW (localization) 

~ 0 


and 


A = 


/~5~ G 5 / 6 Mc /6 
V 96 7 r 2 / 3 c 3 / 2 r 


( 10 ) 


The quantity R is a complicated function of four angular variables (I, ip, 0, (p), and we cannot solve them separately 
only with a single interferometer. Even if the direction (0, <j>) is known ( e.g. from EMW observation), we need at 
least two not-aligned interferometers to solve V, ip and r from observed amplitudes and the relative phases (see 0 
and references therein). 

Now we assume a situation that the three dimensional position (r, 0, (p) of a coalescing binary is determined through 
the afterglow of SHB associated with coalescing (ispiral) gravitational waves detected by LIGO. We start with the 
case for a single interferometer. With gravitational wave observation the chirp mass M c can be determined well by the 
frequency evolution |l(l . Thus we can estimate the intrinsic amplitude A oc M^ 6 /r by combining the chirp mass and 
the distance r from EMW localization. From gravitational wave data we can also get the combination ( AR ) as the 
observed amplitude. We describe the method to constrain the parameter V by dealing with these observed quantities. 
Basic aspects are summarized in table 1. Firstly, we calculate the maximum value of R(V,ip,9,cp) for the observed 
direction (8,<p) of SHB. This is realized with the face-on configuration with \V\ = 1. In this case the polarization 
angle ip does not have meaning due to the symmetry of the geometry, and Rmax is independent on ip, 

R max = R(l,ip,8,<p). (11) 

Secondly, we take the ratio Q of the observed gravitational wave amplitude ( AR ) to the estimated combination 
A • Rmax &S 


Q = 


(AR) 
A ' Rma 


R 


Rr, 


( 12 ) 


The ratio Q is written by 


Q = y -(1 + V 2 ) 2 cos 2 (2 ip + 7 ) + V 2 sin 2 (2 ip + 7 ). (13) 

with 7 = arctan(a 2 /ai). This expression is valid also in the limit (9, </>) —> (7 t/2 ,7 t/4 ) where both the denominator 
and the numerator in eq. m vanish, as the interferometer becomes insensitive to the incident gravitational waves. 
In figure 2 we show the profile of Q with 9 = <p = 0. The ratio Q takes a same value for ±V and, is periodic along 
the '0-direction with period n/2. For a different combination (0, <p) the profile Q in figure 2 shifts to the ^-direction 
as characterized by the angle 7/2 in ea. (11311 . Therefore, without loss of generality, we can study the constraint on 

V for arbitrary (0, (p ) only using figure 2 with which the constrained region for V can be easily read with a given 
Q. Note that the modes higher than the quadrupole order might slightly change this contour map and also affect 
the estimation of the quadrupole amplitude A. These aspects are not peculiar to our method, but common to data 
analyses for gravitational wave astronomy for coalescing binaries. 

For a given value of Q , the maximum of V is realized with the choice ip = 0 in figure 2, and the minimum is with 
ip = n/A. In figure 3 (solid curves) we show the constrained region for V as a function of Q. The upper curve is 

V = Q and the lower one is V = y/2Q — 1. Thus the parameter V is constrained by y / max(2Q — 1, 0) < V < Q. For 
the ratio Q close to 1, the binary must be almost face-on, as this is the only configuration to realize the maximum 
amplitude Q = 1. In contrast, for Q < 1/2, the parameter V can take the range [0, Q\. The allowed region for V is 
very small for Q ~ 1, when the binary is nearly face-on and SHB is expected to be luminous with collimated jet-like 
structures. Even with a smaller Q ( e.g. Q = 0.6) we can correctly discriminate that the SHB is completely off-axis. 

So far we have studied case with single interferometer. By using three widely separated interferometers the direction 
( 0 , (p) of a binary can be determined from the time delays of the gravitational wave signals, and the angles (/, ip) and 
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FIG. 2: The contour plot for Q as a function of V and if> with 6 = 0 and <j> = 0. Contour levels correspond to 0.1 (black) to 
0.9 (white). Q is periodic along r/j-direction with period 7t/2, and Q = 1 for V = 1 (face-on configuration). The contour plot 
Q shifts to ^-direction for different combinations ( 9,cf >). 


the distance r can be also estimated only from gravitational wave observation 0 - As commented earlier, we can, in 
principle, solve two orientation angles (I, ip) and the distance r by using two not-aligned interferometers for a binary 
with known direction. We will see this in §IV. 

LIGO has essentially two interferometers. One is at Hanford, Washington, and another is at Livingston, Louisiana. 
While the separation between them is ~ 3000km (corresponding to ~ 27° measured from the center of the earth), they 
are nearly parallel to increase correlation of burst signals 0 - As a result, it is not easy to observe two polarization 
modes separately and thereby estimate the parameter V = cos I well with using two LIGO interferometers, compared 
with using combination such as LIGO-VIRGO network. In §IV we will return to analyses with multiple interferometers. 
But our studies for a single interferometer would provide rough outlook for determination of the inclination with nearly 
aligned interferometers as LIGO. 


III. OBSERVATIONAL ERROR 


In actual observation, we cannot determine the ratio Q without error A Q. Here we analyze its effects for constraining 
V. As we discussed, the ratio is obtained by three observed values A , R ma x and (AR) (see table 1). We can formally 
write down the relative error for Q as 


A Q A (AR) AA A R max 
Q (AR) A R max 


(14) 


We are dealing with a situation when the afterglow of SHE is observed and the direction ( 9 , <p) is determined very 
well ( e.g . larcsec level). Therefore, the term AR, max /R rnax is negligible compared with other two terms. 

If the orbital precession of the binary is not significant, the estimation error for the observed gravitational wave 
amplitude (AR) has little correlation with errors for other parameters related to the phase, such as, the chirp mass. 
Thus we have 


A (AR) 

TarT 


(SNR)- 1 


(15) 


with SNR being the signal to noise ratio for the detected gravitational wave signal [3- The typical detection 
threshold for a coalescing binary is ~ 8. If we use the temporal information from the observed SHE, the threshold 
somewhat decreases (but less then a factor of 2) {2Qj. 
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FIG. 3: Constraint for V = cos I with given ratio Q. Without error A Q = 0, the parameter V must be within the shaded 
region determined by two solid curves V = Q and V = max(2Q —1,0) . These two curves are obtained by cutting the surface 
in figure 2 along ip = 0 and ip = 7t/4. The dotted curves show the allowed region with 5% error for Q (A Q/Q = 0.05), and the 
dashed curves show the region with 10% error (A Q/Q = 0.10). With a nearly face-on configuration we have AV G [1 — A Q, 1], 


For the intrinsic amplitude A we need the distance r from EMW observation and the chirp mass M c from gravita¬ 
tional wave data. The latter can be determined very accurately, as it is the primary parameter for the time evolution 
of the gravitational wave phase whose information is crucial for detecting gravitational waves with matched filtering 
method 00. In this paper we neglect estimation error for chirp mass. Under this prescription our results do 
not depend on the order of the post-Newtonian expansion in the restricted post-Newtonian approach without orbital 
precession. For a very nearby SHB we might estimate its distance r with various astronomical data. If we use the 
redshift-distance relation (Hubble law; r = cz/Hq at low redshift) to convert its observed redshift z to the distance r, 
the estimated distance r might be significantly contaminated by the peculiar velocity of the host galaxy. At distance 
r ~ 300Mpc (NS+NSs detectable with LIGOII) the uncertainty of the Hubble parameter Hq could be a problem, but 
those of other cosmological parameters (e.g. the density parameter fin) would not be important at these distances. 
The recently reported value Hq = Tl^km/sec/Mpc by WMAP team |2l| contains ~ 5% error. As the distance error 
would be the dominate source of the error A A, we have 


A A A r 
A r 


(16) 


The error for the estimation of the gravitational wave amplitude (AR) and the distance r from EMW observation 
would be independent. Then, using eqs. C3HEHZ3 and putting A R m ax/Rmax = 0, the typical value for relative rms 
error of Q is given as 


A Q 


rms 


Q 



(typical case). 


(17) 


Next we study how the estimation error A Q changes the allowed region for V. For a given Q we expand its possible 
value in the band [max(0, Q — A Q), min(l, Q + AQ)] and solve corresponding V in (V, i/^-plane (figure 2). Then we 
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obtain the following result 


y / max{2 (Q — A Q) — 1, 0} < V < min(l, Q + A Q). (18) 

In figure 3 we added the allowed region for V from the observed value Q whose relative error is A Q/Q = 0.05 (dotted 
curves) and 0.10 (dashed curves). For a almost face-on binary (1 — V <C 1), we have V £ [1 — A Q, 1], With signal to 
noise ratio SNR ~ 10 and the distance error less than 5%, the constraint becomes V £ [0.87,1]. 


IV. MULTIPLE DETECTORS 
A. Solving degeneracy 

In this subsection we discuss observational analysis to estimate the inclination V separately from the polarization 
direction r ip with using networks of detectors. Following Ref.fl2] we introduce two functions er(TV) and e(TV) that 
characterize sensitivity of a network to gravitational waves with various directions and polarizations. As concrete 
examples of networks, we consider the following three combinations; (i) LIGO-Hanford (4km), LIGO-Livingston and 
VIRGO, (ii) LIGO-Hanford (4km) and VIRGO, (iii) LIGO-Hanford (4km) and LIGO-Livingston. The total numbers 
of detectors Nd are three for case (i) and two for cases (ii) and (iii). For simplicity we assume that all detectors have 
identical sensitivity. The total signal to noise ratio for a binary with directional vectors (TV, L), distance r and a fixed 
chirp mass is written as 


SNRl t 


= ^a(TV) [cq(V) + e(TV)ci(V) cos(4$] , 


(19) 


where ip is the polarization angle of the binary angular momentum L measured from a preferred direction that is 
determined by configuration of each network and also fixes two orthogonal polarization bases in the present analysis. 
The combination rg/r represents the signal to noise ratio of single detector for a face-on (|V| = 1) binary at perpen¬ 
dicular direction to the detector plane (6 = 0 in figure 1). For a binary neutron stars with chirp mass 1.2M© the 
distance rg corresponds to ~5000Mpc for advanced LIGO. Here we used numerical results given in |22j based on the 
noise curve for advanced LIGO with its wide band setting {2£j. Using relative sensitivities of LIGO and advanced 
LIGO to double neutron stars ( e.g. [24]), the distance rg for LIGO is ~ 320Mpc. Note that these distances are for 
SNR = 1 with an optimal configuration of a binary and a single interferometer. 

The function cr(TV) shows sensitivity of a network to gravitational wave from direction TV with averaged polarization, 
and generally takes values in 0 < cr(TV) < Nd/2 (Nd'. number of detectors with identical sensitivity). For our three 
networks its maximum value is 1.04 for case (i), 0.65 for case (ii) and 0.94 for case (iii), while its minimum value is 
0.148, 0.118 and 0.0040 respectively. As two LIGO detectors are nearly aligned, they are not effectively complimentary 
to decrease the blind directions with er(TV) -C 1. 

The function e(TV) (0 < e < 1) represents asymmetry of a network with respect to its sensitivity to two orthogonal 
polarization modes coming from direction TV. For a direction TV with e(TV) ~ 1 a network is sensitive only to one 
polarization mode, while two modes are measured with a similar sensitivity for e(TV) ~ 0. With a single detector we 
identically have e(TV) = 1. For our three networks, this parameter can take the maximum value e(TV) = 1 for some 
directions. 

In eq. (HU) two functions Cg(V) and C\(V) are defined by 


co 00 


(1 + 1 / 2)2 
4 


+ V 2 , 


Cl 00 


0 + vY_ V 2 

4 


( 20 ) 


For a face-on binary with |V| = 1, we have ci = 0 and the expression 1191) does not depend on the polarization angle 
ip, as expected. 

Next we evaluate magnitudes of estimation errors A on for fitting parameters at = (r, V, ip) that are related to 
amplitude of the gravitational waveform. We use the Fisher matrix approach that basically uses linear responses of 
the waveform to variations of fitting parameters. In other words, the first derivatives of the waveform by the fitting 
parameters are used to evaluate the expected errors. Under this approach, correlation between the parameters a.i and 
those related to frequency evolution are practically negligible in the restricted post-Newtonian approach 0 

With using eqs.(4.38) and (4.39) in Ref. 0 and integrating out the information of the phase constant ip c (see 
arguments following eq.(2)), the covariance of estimation errors for parameters ( r,V,ip ) are given as 

(AipAip) = E [l + 6U 2 + V 4 + (1 — V 2 ) 2 e(N) cos(4i/S)] /(I — V 2 ) 2 


(21) 
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l^AxpAr) 

= —2Ere(N) sin^^) 

(22) 

(APpAV) 

= 0 

(23) 

(ArAr) 

= 4£r 2 (l + R 2 - (1 - V 2 )e(N) cos(4 ■ip)) 

(24) 

(ArAV) 

= 2ErV(3 + V 2 - (1 - V 2 )e(N) cos(Aip)) 

(25) 

(ARAR) 

= E[ 1 + 6R 2 + R 4 + e(N) cos(4^)(l - R 2 ) 2 ] 

(26) 


r 2 

E ~ r 2 a(N)( 1 - e(7V) 2 )(l - R 2 ) 2 ' 

(27) 


Equation <G3J for the distance error is same as eq.(4.41) in Ref. E2- In the next subsection we will compare the esti¬ 
mation errors for the inclination under various situations, and hereafter, we denote the expression 1261 ) as (AV AV) N 
for notational clarity. The errors become large for binaries close to face-on through the denominator of the factor 
E. This is because the derivatives of the gravitational waveform with respect to the parameters r and R become 
linearly dependent at |R| = 1, and they are degenerated at parameter fitting 0 . In the same manner the error 
becomes large for sky directions with e(N) ~ 1, as we can measure only one polarization mode that are not enough to 
determine parameters (r, V. ill) separately. However, without depending on the values e(N) and R, the combination 
SNRtot oc [(co(R) + e(N)a(N) cos(4^))] 1 / 2 /r can be determined well and its error is given only by the the total 
signal to noise ratio as A (S N R to t) / (S N R to t) = ( SNR to t ) _1 - The situation is same as in the previous sections for 
single detector. We can also confirm this with directly using equations from 1211 ) to ( 1271 . 

Now we study the case that the distance r to a binary is given from EMW observation, assuming its error is 
negligible. After straightforward calculation with eqs. (1211 to (1271 . the estimation errors for parameters {•ip, V) are 
given by a 2 x 2 Fisher matrix, and the V — V component becomes 


(AV AV) r 


r 2 ( 1 — e(N) 2 cos 2 (4 ip)) 

r$a(N)(l - e(N) 2 )(l + V 2 - cos(4^)e(IV)(l - V 2 ))' 


(28) 


Here the suffix r represents the condition that the distance r of the binary is given. Note that the singular behavior 
at \V\ = 1 disappeared. This is because we do not need to solve r and V simultaneously only from gravitational wave 
observation. However, the error (AV AV) X J 2 for the inclination still becomes large with e —> 1. To deal with the 
situation we introduce a new parameter Q defined by 


ARtot = ( cp(V) + e(N)ci(V) cos(40) \ 1/2 = / c 0 (V) + e(N) Cl (V) cos(4VQ \ 1/2 
ARmax V c o(l) + e(iV)ci(1) cos(4'0) ) V 2 / 


(29) 


and change variables from (•ip, V) to ( Q , V). The parameter Q for multiple detectors is a simple generalization of the 
previous parameter Q in eq.(12) defined for single detector. In figure 4 we show the allowed region of the combination 
( Q , V) for given network parameter e. For a direction with e = 1 the allowed region coincides with that for Q-V given 
by the shaded region in figure 3. When we decrease e form 1 to 0, the region shrinks as in figure 4, and it becomes a 
single curve Q = (co(V)/2) 1 / 2 at e = 0. 

Next we consider a Gaussian probability distribution function (PDF) P(Q X , V x \Q, V) of the estimated parameters 
(Qx, V x ) around their true values ( Q, V). We calculate the covariances of errors for parameters ( Q, V) by using the 
2x2 Fisher matrix for the original combination (-0, V). Even if the error (AVAV) r ' is large at e(N) ~ 1 as shown 
in eq. |28|) . the parameter Q can be estimated relatively well, as in the case for single detector. For example, with 
a 1 — a error ellipsoid that is stretched toward R-direction as shown in figure 4, we can estimate the inclination V 
much better by limiting the combination (Q x , V x ) into its allowed region discussed above. This difference shows one 
of potential problems of a simple evaluation based on Fisher matrix without dealing with global domain of the fitting 
parameters. Therefore, we define a new PDF Pc(Q x , Vx\Q, V) with the following two steps. We firstly narrow the 
parameter space ( Q X ,V X ) of the original Gaussian distribution P(Q X , V X \Q, V) using the Q-V constraint in figure 4, 
and secondly renormalize the PDF appropriately. Then we calculate the estimation error (AVAV) c as 


(AVAV) C = J dV x dQ x P c (Q x , V X \Q,V)(V X - V) 2 . 


(30) 


Here the suffix C represents to use the constraint. 
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FIG. 4: Constraints for V and Q for given parameter e(N) that represents asymmetry of the sensitivity of a network to two 
polarization modes coming from direction N. For e = 1 the combination (Q, V ) is in the region bounded by two solid curves. 
This region is same as the shaded region in figure 3 for single detector. The short-dashed curves are for the boundary of the 
allowed region t = 0.5. For t = 0 the allowed region is on the long-dashed line Q = \cq(V) /2] 1 ' 2 . The long ellipsoid around 
(Q,V) = (0.9, 0.9) is an example of 1 — a error ellipsoid described in the main text. 


B. Monte Carlo analysis 


In this subsection we compare estimations errors AV given by (AV AV) 1 ^ 2 , (AVAV)l /2 and ( AVAV)^ 2 presented 
in the previous subsection. The first one (AV AV) ^ is the prediction by Fisher matrix for a binary whose distance 
is not given from EMW observation. The second one (AV AV) l J 2 is also predicted by Fisher matrix but for binaries 
with known distances. For the third one we apply the Q-V constraint to binaries with known distances. For each 
network (i)-(iii) we prepare sample of “detectable binaries” with SNR to t > 8 in the following manner. From eq. m 
the maximum distance r max of binaries with the threshold SNR to t = 8 is given as r max = ro[2max<j(iV)] 1 / 2 /8. A 
binary with SNR tot > 8 must be in a sphere with radius r max around the Earth, since this distance r max is given 
for optimal configuration. In this sphere we put a binary with random position and orientation. If its SNR to t is 
larger than 8, we add it to our sample of “detectable binaries”. We continue this until the total number of our sample 
becomes 200 for each network (i)-(iii). 

We find that, out of these samples, the numbers of binaries with \V\ > 0.8 are 102 for case (i), 99 for case (ii) and 
101 for case (iii). For binaries with \V\ > 0.9 we have 64, 62 and 60 respectively. While the intrinsic distribution of 
the parameter \V\ is homogeneous in the range [0,1] with random direction N and orientation L , our “detectable 
binaries” have skewed distributions of \V\ after selected with a signal-to-noise threshold. 

In figure 5 we show the distribution of the estimation errors AV for case (i). The upper panel shows the results for 
all 200 binaries and the lower panel is only for 102 binaries with \V\ > 0.8. Similar results are given in figs. 6 and 

I /n 

7 for cases (ii) and (iii). Comparing distributions for two Fisher matrix predictions (AV AV) r (dashed curves) and 

(AVAV)x (dotted curves), it is apparent that the estimation of the parameter V is significantly improved, if the 
distance r to the binary is determined by EMW observation. The differences between two curves become larger for 
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FIG. 5: cumulative distribution of the estimation errors AV for binaries with using two LIGO and VIRGO detectors (case 
(i)). The solid curves are results for (AVAV)J/ 2 (known distance from EMW observation and Q-V constraint). The dashed 
curves are given by the Fisher matrix approach (AVA V)\! 2 (known distance from EMW observation). The dotted curves are 
for (AVAV)J/ 2 (unknown distance from EMW observation). The lower panel is for binaries with \V\ > 0.8. 


subsample with \V\ > 0.8 due to the degeneracy for the parameter fitting at the face on limit |V| —> 1. 

We can study the effects of the Q-V constraint with solid curves and dashed curves in figs.5-7. Note that this 
constraint does not always decrease the estimation error AV, as shown in figs. 5 and 6. However, the tail of large 
AV are removed by the constraint. As a result, the error AV becomes smaller than ~ 0.25 for our samples. The 
tails are mainly made by binaries with e ~ 1 for which it is difficult to observe two independent polarization modes 
and thereby solve two parameters (V, ?/>) separately. The constraint works well especially for case (iii) with two LIGO 
detectors whose orientations are nearly aligned. In this case the situation is similar to the observation with using a 
single detector as studied in §2. In the lower panel of figure 7 the median value of the error AV is reduced by ~ 40% 
by using the constraint. This indicate that we should be careful to use a simple Fisher matrix approach for evaluating 
the parameter estimation error for inclination. 


V. DISCUSSIONS 

We discussed prospects of a method to constrain the inclination angle of a coalescing compact binary by detecting its 
gravitational waves associated with a three-dimensionally localized SHB. With this method we can get an important 
geometrical information to understand the properties of SHBs and their afterglows as a function of the viewing angles 
of the jets (see e.g. j2f| for long-soft bursts). We should comment that there would be a selection effect toward larger 
\V\ for a simultaneous detection of gravitational waves and SHB. This is because the observed gravitational wave 
amplitude would be larger (see §IV.B) and the SHB would be also luminous with the face-on configuration Iff, I20I. 

By analyzing nearby samples with our method we might establish an efficient empirical criteria to select almost 
face-on binaries using observed properties of EMW signals. Then we can study cosmological parameters with face- 
on binaries at relatively high redshift where cosmological effects beyond the Hubble-law become important. The 
following is the outline of this approach. We first obtain the intrinsic amplitudes A by putting R = R m ax for the 
observed gravitational wave amplitude (AR) for binaries that are likely to be nearly face-on according to the empirical 
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FIG. 6: Same as figure 4 with LIGO-Hanford and VIRGO detectors (case (ii)). 




FIG. 7: Same as figure 4 with two LIGO detectors (case (iii)). 
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criteria. Then the luminosity distances r are obtained from eq.(10) and the estimated chirp masses. In this manner 
we might observationally study the redshift-distance relation with the redshift information naturally obtained from 
identified host galaxies with EMW observation. Note that this argument to estimate the distance r is different from 
the previous method in this paper to constrain the inclination angle by simply converting the observed redshift to the 
distance r for nearby binaries. Considering the strong correlation between the inclination V and the distance r for 
parameter estimation 0 , this could be a powerful approach to investi gat e the dark energy with future gravitational 
wave detectors that can detect compact binaries at z ~ 1 (see also [T6| ]2£ [2(|). 
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